Decay of energy and suppression of Fermi acceleration in a dissipative driven 

stadium-like billiard 
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The behavior of the average energy for an ensemble of non-interacting particles is studied using 
scaling arguments in a dissipative time-dependent stadium-like billiard. The dynamics of the system 
is described by a four dimensional nonlinear mapping. The dissipation is introduced via inelastic 
collisions between the particles and the moving boundary. For different combinations of initial 
velocities and damping coefficients, the long time dynamics of the particles leads them to reach 
different states of final energy and to visit different attractors, which change as the dissipation is 
varied. The decay of the average energy of the particles, which is observed for a large range of 
restitution coefficients and different initial velocities, is described using scaling arguments. Since 
this system exhibits unlimited energy growth in the absence of dissipation, our results for the 
dissipative case give support to the principle that Fermi acceleration seem not to be a structurally 
stable phenomenon. 

PACS numbers: 05.45. Pq, 05.45.Tp 



Some dynamical properties of a dissipative 
time-dependent stadium-like billiard are stud- 
ied. The system is described in terms of a four- 
dimensional nonlinear mapping. Dissipation is 
introduced via inelastic collisions of the particle 
with the boundary, thus implying that the parti- 
cle has a fractional loss of energy upon collision. 
The dissipation causes substantial modifications 
in the dynamics of the particle as well as in the 
phase space of the non-dissipative system. In par- 
ticular, inelastic collisions are an efficient mecha- 
nism to suppress Fermi acceleration of the parti- 
cle. We show that a slight modification of the in- 
tensity of the damping coefficient yields a change 
of the final average velocity of the ensemble of 
particles. Such a difference in the final plateaus 
of average velocity is explained by a large num- 
ber of attractors created in the phase space by 
the introduction of dissipation in the system. We 
also described the behavior of decay of energy 
via a scaling formalism using as variables: (i) ini- 
tial velocity (Vb); (ii) the damping coefficient (7) 
and; (iii) number of collisions with the boundary 
(n). The decay of energy, leading the dynamics 
to converges to different plateaus in the low en- 
ergy regime, is a confirmation that inelastic col- 
lisions do indeed suppress Fermi acceleration in 
two-dimensional time-dependent billiards. 



I. INTRODUCTION 

Billiard problems symbolize the dynamics of a point- 
like particle (the billiard ball) which moves inside a com- 
pact region Q, that in the context of the mathematics 



of billiards is known as a billiard table (or for short, bil- 
liard). Inside the billiard, the particle moves in straight 
lines until it reaches the boundary where the specular re- 
flection rule is used (i.e., mirror like). This implies that 
the incidence angle is equal to the reflection angle upon 
collision: for a general reference, see, e.g., [lj]. The first 
studies about billiards began with Birkhoff in 1929 Q, 
who proposed the billiard ball motion in a manifold with 
an edge. After the pioneering results of Sinai [|J , Buni- 
movich 0, [H and Gallavotti @ , who gave mathematical 
support to the field, several applications of billiards have 
been found in different areas of research including optics 
-Tl0|| . quantum do ts [l~Tl| . microwaves [Hj], astronomy 
laser dynamics [14j and many others. 

In the case when the billiard boundary is time- 
dependent, dQ = dQ(t), depending on the shape of the 
border, the particle can accumulate energy under the ef- 
fect of successive collisions, leading to the phenomenon 
known as Fermi acceleration (FA). Introduced in 1949 
by Enrico Fermi (Tol l as an attempt to explain the high 
energy of the cosmic rays, FA basically consists in the 
unlimited energy growth of a point-like particle suffering 
collisions with an infinitely heavy and time-dependent 
boundary. One of the important questions which arises 
from the study of 2-D time-dependent billiards is: what is 
the condition which leads the particle to experience unlim- 
ited energy growth? As previously discussed in Ref. [l6| . 
the Loskutov-Ryabov-Akinshin (LRA) conjecture claims 
that the chaotic dynamics for a particle in a billiard with 
static boundary is a sufficient condition to produce FA 
if a time-perturbation to the boundary is introduced. 
Later, F. Lenz et al in Ref. [I?} studied the case of a 
specific time-perturbation to the elliptical billiard, which 
is integrable for the static boundary. In this elliptic par- 
ticular case, FA is produced by orbits that "cross" in the 
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phase space the region of the separatrix, which marks the 
separation from motion of two kinds: (i) libration and 
(ii) rotation; therefore characterizing a change in the dy- 
namics from librator to rotator 

0113 > 

or vice- versa. This 

"crossing" behavior is indeed assumed to be the FA pro- 
duction mechanism. Latter on [Lsl ] was shown that when 
the crossings were stopped, the FA is suppressed and the 
energy of the particle is constant for long time dynamics. 
As illustrated in pj} and confirmed in [l8[ for a differ- 
ent time-perturbation, the existence of a separatrix curve 
in the phase space, which is observed in the static case, 
turns into a stochastic layer when a time perturbation 
is introduced on the boundary, and produces the needed 
condition for the particle to accumulate energy along its 
orbit leading the dynamics to exhibit FA. 

While the condition to produce FA in billiards is well 
understood, a question which naturally arises is: "What 
should one do to suppress FA in time-dependent billiards 
if the energy of the particle is unlimited?" The study of 
suppression of FA is quite recent for two-dimensional bil- 
liards and one can consider different kinds of dissipative 
forces. To illustrate a few of them, the introduction of 
inelastic collisions has been discussed for the oval-like bil- 
liard [HI, Lorentz gas [H and elliptical billiard QJ]. The 
introduction of a drag force can also be considered a way 
to suppress FA, as demonstrated in the oval-like billiard 
[U and elliptical billiard [H|. 

In this paper, we consider the dynamics of a time- 
dependent stadium-like billiard, aiming to understand 
and describe the behavior of the average velocity for an 
ensemble of initial conditions when dissipation is intro- 
duced. First, we construct the equations describing the 
dynamics of the model, including those governing the 
reflection rule. Then we investigate the dynamics for 
regimes of high and low energy. The initial high energy 
regime, as will be shown, decays in time according to a 
power law. Critical exponents are derived and scaling 
arguments are used to describe a scaling invariance for 
the average velocity in the regime of high energy. Gener- 
ally the dynamics of the system depends on the control 
parameters including those controlling the non-linearity 
of the system. As they are changed, average quantities 
of some observables exhibit typical behavior observed in 
phase transitions [23|. Near the phase transition, crit- 
ical exponents can be defined and a scaling investiga- 
tion can be carried out. Our results for the dissipative 
stadium-like billiard show that the phenomenon of FA 
is suppressed even for small dissipative coefficients, then 
changing the regime of unlimited energy growth to lim- 
ited growth. Since the conservative version of this bil- 
liard presents FA [23 - |26j . this result give support to the 
principle that FA seem not to be a structurally stable 
phenomenon 27j. 

This paper is organized as follows: In Sec. [II] we con- 
struct the equations that describe the dynamics of a dis- 
sipative time-dependent stadium-like billiard. Section lTOl 
is devoted to discuss our results and is therefore divided 
in two parts: (i) The first one includes the investigation 



of the chaotic transient considering high and low energy 
regimes, including the scaling investigation and critical 
exponents, (ii) In the second part, the final convergence 
for the velocity is studied as function of the dissipation 
and initial energy regime. Our final remarks and conclu- 
sions are drawn in Sec. IIVI 



II. THE MODEL AND THE MAPPING 

In this section we construct the equations that describe 
the dynamics of the system. The model describes the 
dynamics of a point-like particle suffering inelastic colli- 
sions with a time-dependent stadium billiard. Inelastic 
collisions are introduced by two distinct damping coef- 
ficients 7 £ [0,1] and (3 £ [0,1], where 7 corresponds 
to the restitution coefficient with respect to the normal 
component of the boundary at the instant of the colli- 
sion while (3 is the restitution coefficient with respect to 
the tangential component. For 7 = j3 = 1, as expected, 
results of the non-dissipative case are all obtained. To 
construct the dynamics of the model, we have to consider 
two distinct situations: (i) successive collisions and; (ii) 
indirect collisions. For case (i), the particle suffers suc- 
cessive collisions with the same focusing component. On 
the other hand, in (ii), after suffering a collision with a 
focusing boundary, the next collision of the particle is 
with the opposite one where the particle can, in princi- 
ple, collide many times with the parallel borders. We 
have considered that the time dependence in the bound- 
ary is R(t) = Rq + r s'm(wt) , where Ro ^> r. The velocity 
of the boundary is obtained by 



R{t) = B(t) = B cos(wt) 



(1) 



with Bo = rw and r is the amplitude of oscillation of the 
moving boundary while w is the frequency of oscillation. 
In our approach, the dynamics of the model is described 
in its simplified version in the sense that the boundary is 
assumed to be fixed, but, at the moment of the collisions, 
it exchanges energy with the particle as if it was moving. 
We consider fixed w = 1. 

The dynamics is evolved considering the variables 
(a n , (fin, t n , V n ), where a is the angle between the trajec- 
tory of the particle and the normal line at the collision 
point, tp is the angle between the normal line at the col- 
lision point and the vertical line in the symmetry axis. 
We assume that V is the velocity of the particle and t is 
the time at the instant of the impact while the index n 
denotes the n th collision of the particle with the bound- 
ary. As an initial condition, we assume that at the initial 
time t = t n , the particle is at the focusing boundary 
and the velocity vector directs towards to the billiard ta- 
ble. In the notation, all variables with (*) are measured 
immediately before the collision. Figure [1] shows an il- 
lustration of a trajectory with successive collisions, where 
R = (a 2 + 4b 2 ) /(8b) and $ = arcsin(a/2i?). The control 
parameters a, b and I are drawn in Fig|2] If a = 2b the 
original stadium billiard is recovered. 
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FIG. 1. (color online) Dynamical variables and a typical tra- 
jectory with successive collisions. 

Let us consider case (i) first. For the occurrence of 
a successive collision, it is necessary that | <p n +i |< 
Therefore, according to Fig. [T] and the specular reflec- 
tion, we have that 

* — r\ 

a n+l — a n i 

ifn+i = <pn + 7T - 2a n (mod 2-k) , 
_ 2iT , cos(q ra ) 

Considering the case of indirect collisions, case (ii), it is 
necessary that | (p n +i > To obtain the equations de- 
scribing the dynamics, it is useful to consider the unfold- 
ing method Then two auxiliary variables, ip> which 
is the angle between the trajectory and the vertical line 




FIG. 2. (color online) Illustration of the: (i) unfolding 
method; (ii) control parameters and (iii) a typical trajectory 
with indirect collisions with the focusing component. 



at the collision point and, x n , which is the projection of 
the particle position under the horizontal axis are used. 
From geometrical considerations of Fig. [T]we obtain that 
•ip n = a n — tp n . One sees also that x n is the summation of 
the line segments AB + BC + CD. Taking into account 
the expression of ip n , and after some algebra, we obtain 
x„ = cos (^„) [sin(a n ) + sin(<I> - ijj n )\. The recurrence re- 
lation between x n and x n +\ is given by the unfolding 
method, described in Fig. [3J as x n+ i = x n + I ta,n(ip n ) . 

To obtain the angular equations, we invert the par- 
ticle motion, i.e., consider the reverse direction of the 
billiard particle, then the expression that furnishes x n is 
also inverted and the angle a„ becomes a* . Resolving it 
with respect to a* , taking into account that this angle is 
changed in the opposite direction, then the angles a n and 
ip n must have their sign reversed. Moreover, the incident 
angle a* assumes when the motion of the particle 

is re-inverted. The expressions of </?„+i and the time t n+ \ 
are obtained by simple geometrical considerations of Fig. 
[2] Thus, we obtain the mapping for the case of indirect 
collisions as 

a n+i = arcsin[sin('0 n + $) - x n +i cos(ip n )/R] , 

fn+l = - Ol* n+1 , 

i?[cos((/? ra ) + cos(y?„ + i) - 2cos($)] + I 

V n COS(ipn) 

For both cases (i) and (ii) , the recurrence relations for 
V n and a n are the same. Let us discuss how to obtain 
them. Expressing the two components of the velocity 
vector before the collisions we have that 

V n ■ T n+ i = -Vn COs(a* ) , 

V n ■ N n +x = V n sm(a* n ) . 

Given that the moving boundary is not an inertial refer- 
ential frame, we assume that at the instant of the colli- 
sion, the wall is instantaneously at rest, then the reflec- 
tion laws are given by 

V n+ i ■ T n+ i = /3V n ■ T n+ i , 

VUl ■ N n+ i = ■ N n+ i , 

where /? and 7 are respectively the restitution coefficients 
with respect to the tangent and the normal components 
of the motion. We stress that V^ +1 -T n+ i and V^ +1 - N n +i 
are the components of the velocity of the particle mea- 
sured in the referential frame of the moving wall with 
respect to the tangent and the normal components re- 
spectively. 

In the inertial referential frame, we have that the equa- 
tions for the components of the velocity are given by 

V n +i ■ T n+ i = f3V n ■ T n+ i + (I — j3)B(t n+ i) ■ T n+ i , 

V n+ i ■ N n+1 = -7% ■ Nn+1 + (1 + l)B{t n+1 ) ■ N n+1 , 

where B(t n +i) is obtained from Eq.([T]) evaluated at the 
time t n+ i. Finally, the expression of V n +\ is given by 

I V n +i |= \j (K+i • T„+i) 2 + (V n+1 ■ AT„ +1 )2 . (2) 
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The last equation refers to the reflection angle, which 
according to the reflection law is given by 



Vn 



*n+l 



arcsm 



Vn 



■sin« +1 ) 



(3) 



+1 



III. 



DISCUSSIONS OF THE DYNAMICS AND 
NUMERICAL RESULTS 



To investigate the dynamics of the model, we set as 
fixed the parameters a = 0.5, b = 0.01, I = 1 and the 
amplitude of oscillation of the moving wall as Bq = 0.01. 
The reason for keeping the parameters fixed is because 
some observables are scaling invariant with respect to the 
control parameters, as discussed for the static case in Ref. 
[28| . Therefore for this paper we chose to investigate the 
dynamics by the variation of the parameter 7 as well as 
the initial velocity Vq. The parameter f3 is considered 
fixed as f3 = 1. Moreover we stress that collisions of 
the particle with the straight segments of the border are 
considered elastic. 

Figure [3] illustrates the dynamics of the particle on 
the variables (tp,£), where £ = x/a, for the conserva- 
tive case and the dissipative dynamics respectively. It is 
known from the literature that, when dissipation is intro- 
duced in the dynamics, the invariant curves surrounding 
the fixed points might be destroyed and the elliptic fixed 
points turn into sinks [29] ], Specific discussions of fixed 
points in the non dissipative case can be seen for the 
stadium-like billiard in Refs. [24T - I261 . |28| . The procedure 
used to construct the phase portraits was to consider 25 
different initial conditions and evolve each initial condi- 
tion until 10 7 collisions. For the conservative case, FA 
is indeed oberserved for the stadium billiard for an ini- 
tial velocity higher than the critical resonant one [zij . 
and the phase space for this unlimited energy growth is 
shown in Figj3ja). On the other hand, the darker re- 
gions of Fig. M.b) mark the convergence to the attrac- 
tors while the spread points along the plot identify the 
seemingly chaotic transient. The control parameters and 
initial conditions used in the construction of the figure 
were: (a) V = 10, 7 = 1.0 and (b) V = 10, 7 = 0.999. 
Comparing FigJHJa) and (b), one can see the "birth" of a 
period-2 atractor, and two period one atracting ponits. 

In order to have a better understanding of the dynam- 
ics, particularly the convergence of the initial conditions 
to the attractors, we concentrate to study the dynamics 
and hence some properties of the average velocity of the 
particle. We consider dependence of the average veloc- 
ity as function of both: (a) number of collisions with the 
boundary n; (b) initial velocity Vq and; (c) restitution 
coefficient 7. The average velocity is therefore defined as 



_ 1 M 



(4) 




FIG. 3. Phase portrait of (V>,£) for: (a) Vo = 10, 7 = 1.0 and 
(b) Vo = 10, 7 = 0.999. 



tions (a, tp), from which £ and ip are described, and Vi is 
expressed by 



1 - 

V t (n,V 0ll ) = -y2V k 

71 ' 



(5) 



k=l 



where M is an ensemble of 5000 different initial condi- 



where n is the number of collisions with the moving wall. 

The numerical investigations were carried out in two 
different ways. To understand the behavior of the average 
velocity and hence the energy of the particle at the range 
of large initial velocity we consider two regimes of time: 
(i) short time, mainly marked by the dynamics evolving 
through a transient and; (ii) long time where the dynam- 
ics has already reached the attractor. For high initial 
velocity, the particle experiences a decay in the average 
velocity marked by a transient in the dynamics. More- 
over the decay of energy is described by a homogeneous 
function with critical exponents, (ii) For long time, a 
statistics is made to the regime of convergence ir order 
to understand the role of the dissipation and of the ini- 
tial energy regime. The simulations were evolved up to 
n = 5 x 10 8 collisions and considering an ensemble of 
5000 different initial configurations uniformly chosen as 
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tp £ [0, $] and a £ [0, 7r/2]. Due to the axial symmetry of 
the stadium-like billiard, a negative range of the initial 
conditions is not needed. 



A. Transient for short time 

In this section we concentrate to investigate the initial 
transient. Therefore we consider two ranges of initial ve- 
locity: (i) large and (ii) low. Let us start with the high 
energy. Figures E(a,b) illustrate the behavior of V as 
function of the number of collisions. In Fig. Eta) we as- 
sume as fixed the initial velocity as Vq = 100, and varied 
the restitution coefficient 7. In Figure HJb) we consid- 
ered fixed 7 = 0.999 and varied the initial velocity Vq. 
The range of Vq was chosen in such a way to configure 
a very large initial velocity as compared to the maxi- 
mum boundary velocity Vq » ru. The curves shown 
in Figs. HJa) and Sib) indeed exhibit similar behavior 
for short time. They begin in a constant regime for each 
initial velocity and suddenly, depending on the value of 
the damping coefficient 7, they experience a crossover 
n x marking a change from a constant regime and bend 
towards a regime of decay according a power law. A care- 
ful fitting in the curves gives that the decay exponent is 
(w -1. After the decay observed for large n, the curves 
of V saturate in different plateaus of low energy, which 
may depend of both Vq and 7. The plateaus character- 
ize indeed different attractors to where the dynamics has 
converged to. For large enough time the convergence re- 
gions are around the range V £ (0.07, 0.6) The investiga- 
tion of these plateaus will be made in the next subsection. 

Given that the initial behavior of V for the range of 
large initial Vq is similar even for different values of Vo 
and 7, we can suppose that: 

• (i) V oc V a , for n <C n x , where a is a critical expo- 
nent and n x is the characteristic crossover collision; 



• (ii) V oc (p^) ! for n ^> n x , where £ 
power law decaying exponent; 



-1 is the 



("i) {%) K 
dynamical exponents. 



On the scaling hypothesis (iii) we considered (1 — 7) 
instead of 7, because we want to consider the transition 
(1 — 7) — > + . We see also that different initial velocities 
produce different curves of V but with the same negative 
slope. Therefore the transformation n — »■ ti/Vq makes all 
the curves coincide in the decay regime. This transfor- 
mation together with some curves of V are shown in Fig. 
[5j After this transformation, all curves of V start in a 
constant regime and then decay together as a power law 
with exponent £ w — 1. 

Given that the initial behavior of V is constant, we 
conclude that a — 1. The critical exponents z\ and z 2 
are obtained respectively by power law fits on the plots 
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7) Z2 , where Z\ and z 2 are 



FIG. 4. (color online) Plot of V vs n for the large initial 
velocity. The parameters and initial conditions used were: 
(a) Vb = 100 and different restitution coefficients 7 and; (b) 
different initial velocities and a fixed 7 = 0.999. 



(n x /V ) x V and (n x /V ) X (1 - 7). Figure [S] gives us 
that = -0.99(1) and z 2 = -0.968(1). 

The three scaling hypotheses allow us to describe the 
behavior of V for short n (before the convergence to the 
constant plateau) formally as a scaling function of the 
type 



V(n/V ,V (l-7)) = Ay(A Ql n/Vb,A fol Vb(l-7)) , (6) 

where A is a scaling factor and a± and b\ are scaling 
exponents. Assuming that X ai n/Vo constant, we have 



(7) 



Substituting Eq.([7]) in (|6]), we obtain 

y(n/Vb,Vb(l-7)) = Vo~*Vl(l,A-^Vo(l-7)) . (8) 

If we compare Eq. (|5|) with the hypothesis (ii) and assum- 
ing that V\ is constant for we obtain 



c 



1 



(9) 
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FIG. 5. (color online) Plot of V as a function of: (a) n and 
(b) n/Vo- The parameter 7 was fixed as 7 = 0.999 and 5 
different initial conditions were used, as labeled in the figure. 



and given the constancy of the initial regime, we conclude 
that a = 1, yielding b\ = —1. 

Comparing now Eqs.® with ([TO)) and after straight- 
forward algebra, we obtain that 



(13) 



When Eq. (fT3"|) is compared with the scaling hypothesis 
(iii) we conclude that 



Ol 



a. 

c 



zi = z 2 



(14) 



This procedure and the critical exponents let us rescale 
properly both axis of the V vs n plot and obtain a single 
and universal curve for the short time transient dynam- 
ics, as shown in Fig. [7] Basically, this result confirms 
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FIG. 6. (color online) Plot of: (a) n x /Vo vs Vo for 7 = 0.999 
and; (b) n x x (1 — 7) for Vo = 100. After fitting the data we 
obtain Zi = -0.99(1) and z 2 = -0.968(1). 



and given that the critical exponent ( w — 1, obtained 
by fitting a power law to the decay regime, we have that 
ai = 1. 

Choosing now A bl Vq(1 — 7) constant, we have 



A = (Vb(l- 7 )r~ 
Replacing Eq. dTUl) in ([5]), we obtain 



(10) 



V(n/V , Vo(l - 7)) = V Q (1 - 7rW 2 (A--n/Vo, 1) • 

(11) 

A comparison of Eq. (ITT1) with hypothesis (i), and as- 
suming V% constant for n <C leads to 
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(12) 



(n/V n )/V n '(1-7) 



FIG. 7. (color online) Plot of: (a) V vs n, and (b) overlap 
of the initial transient of all curves of (a) onto a single plot, 
after a suitable rescale of the axis. 

that, independent of the initial velocity and the control 
parameter 7 = 1, the behavior of the transient for the 
high energy regime of V curves is scaling invariant with 
respect to Vo and 7. 

Let us now consider the case where the initial velocity 
is low, therefore the behavior of the transient is different. 
The curves of V exhibit a regime of growth until they 
reach the convergence regions, as shows Figs. [8]Ja,b). 
The combination of control parameters Vo was set in or- 
der to keep the initial velocity in the low energy regime. 
The dissipation and the initial velocity are labeled in the 
figure. 

The presence of the attractors indeed define the region 
to where the curves of V converge to and therefore satu- 
rate. Although one may think this is quite a paradoxical 
behavior in the sense that dissipation leads to a regime of 
growth, if one looks deeper at the dynamics, indeed the 
particle is only converging to an attractor wich is located 
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FIG. 8. (color online) Plot of V vs n for the low energy regime: 
In (a) Vo = 0.05 was fixed and the dissipation parameter 7 
was ranged while; (b) 7 = 0.999 was fixed and the initial 
velocity for the a initial energy regime was ranged. 



FIG. 9. Plot of V x t mod(2n) after the chaotic transient 
for: (a) V = 1000, 7 = 0.99; (b) Vo = 100, 7 = 0.9999; (c) 
Vo = 100, 7 = 0.999 and (d) V = 10, 7 = 0.999. 



at higer energy compared to the initial velocity. Consid- 
ering that this attractor is not at infinite velocity, the FA 
is suppressed. 



B. Sinks, attractors and Convergence of V 

Once the behavior of the chaotic transient for both 
high and low energy regimes is described, let us concen- 
trate our efforts to investigate the convergence regions, 
the role of the sinks and their dependence according the 
dissipation and initial velocities. 

As shown in Fig. [31 the fixed points in the space (i/j, £) 
become sinks. However a visualization of the phase por- 
traits in such a variables does not give any conclusions 
regarding the final velocity. Then it is natural to look at 
the plots of V vs t mod(2n) where t is the time shown 
in Fig. |9l We can see that the convergence regions for 
the curves of V illustrated in Figs. 0] and [5] are around 
V £ (0.07, 0.6). It is also possible to see from Fig. [S]many 
different sets of attracting fixed points and more complex 



attractors. A zoom-in window shows better some of the 
sinks in Figs. |HIa,c). In particular, one can enumerate in 
Fig. [Hid) at least 13 different attractors. Each attractor 
of this set produces a different plateau in the asymptotic 
curves of V. This multitude of attracting points for the 
convergence zone is the reason why a scaling treatment 
for long time is indeed a real challenge. To construct the 
figures, we set 100 different initial conditions, each one 
of them evolved in time until 5 x 10 7 collisions. 

Each attractor has its own influence over the dynamics 
which is dependent on the size of the basin of attraction. 
A way to see this is constructing a histogram of frequency 
of initial conditions showing convergence to a particular 
attractor. Figure [TU] shows the corresponding histogram 
of frequency of visited attractors for the same control 
parameter used in Fig. [5] however with a large ensemble 
of initial conditions, indeed 25 x 10 4 of them where each 
one of them was evolved until 5 x 10 s collisions. 

Figure [TOT a) shows that the most visited attractors 
are those located around V £ (0.08,0.09). A compar- 
ison of this result with Fig. H^a) shows one main at- 
tractor indicating a possibly period-2 attractor. Fig- 
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FIG. 10. (color online) Histogram of frequency for the con- 
vergence velocity region for: (a) Vq = 1000, 7 = 0.99; (b) 
V = 100, 7 = 0.9999; (c) V = 100, 7 = 0.999 and (d) 
Vo = 10, 7 = 0.999. 



lire ITUf b) shows that the distribution of final velocity 
of the orbits, for the combination of control parameters 
7 = 0.9999 and Vo = 100, are more concentrated in 
a range V £ (0.3,0.4). Again, if the results are com- 
pared with Fig. HJb), one sees that the visited region 
V £ (0.3,0.4) corresponds to the seemingly cyclic at- 
tractor. Also, Figs. ITo7 c.d) show many different visited 
regions corresponding to the regions of final velocity as 
shown in Figs. |UJc,d). 

The influence of the attractors arc dependent on the 
control parameters used. Therefore specific attractors 
can be more influence than others for a combination of 
control parameters. In order to understand the depen- 
dence of the attractors with the dissipation parameter 
and the initial velocity, we constructed a histogram of 
frequency taking into account different control parame- 
ters and initial velocities, as shown in Fig. [TTJ After 
a careful look at Fig. [Til we conclude that the attrac- 
tors for "more" dissipative systems, for example the ones 
for 7 e [0.95, 0.999], prefers the regions of lower velocities 
even considering the high and low initial velocity. On the 



other hand the "less" dissipative dynamics, for example 
for the range 7 6 [0.9995,0.999999], prefers regions of 
higher velocities as compared to the previous case. Each 
histogram shown in Fig. [TTJ was constructed consider- 
ing 25 x 10 4 different initial conditions where each one 
of them was evolved up to 5 x 10 8 collisions. Figures 
lllf a.b) represent the histograms for the initial high ve- 
locity while Figs. [TTT c.d) show the histograms for the low 
initial velocity. In particular, one can notice that in Fig. 
[TlTd). the legend used for each parameter of dissipation 
7, applies to all plots of Fig. [TO 
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FIG. 11. (color online) Histogram of frequency for different 
parameters 7. The control parameters and initial condition 
used were: (a) Vo = 100, (b) Vo = 10, (c) Vo = 0.1 and (d) 
V = 0.01. 



IV. CONCLUSIONS 

In summary, we considered a time-dependent stadium- 
like billiard with dissipation introduced via inelastic col- 
lisions. With the introduction of the dissipation, we have 
shown that FA is suppressed. High initial velocities, af- 
ter a crossover time, decay as a power law of the num- 
ber of collisions with the border, while low initial energy 
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regimes, lead the particle to present a regime of "grow- 
ing" to the same convergence region of the orbits with a 
high energy chaotic transient. This chaotic transient for 
the high energy regimes, was characterized as function of 
both Vo and 7. Scaling arguments were used to overlap 
the behavior of V for short n, showing that the dynamics 
of short time is scaling invariant with respect to Vo and 7, 
if we considered high initial energy regime. The system is 
shown to have many attractors, where some of them are 
sinks and others are more complicate. When the damping 
coefficient is varied there is a tendency for more dissipa- 
tion bringing the dynamics to visit more often the region 
Vfi n ai ~ 0.30. On the other hand, less dissipative dy- 
namics prefers the attractors around Vfi na i £ (0.3,0.6). 
Finally, it is clear that the unlimited energy growth is 
interrupted with the presence of inelastic collisions there- 



fore leading to one more example where Fermi accelera- 
tion seems not to be a structurally stable phenomenon. 
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